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C- ' Abstract. The Efroimsky perturbation scheme for consistent treatment of 

^ ' gravitational waves and their influence on the background is summarized and 

compared with classical Isaacson's high-frequency approach. We demonstrate that 
the Efroimsky method in its present form is not compatible with the Isaacson limit 
of high-frequency gravitational waves, and we propose its natural generalization 
to resolve this drawback. 
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1. Introduction 

Some time ago Efroimsky introduced and developed new formalism for a consistent 
treatment of weak gravitational waves [1,2]. This interesting mathematical framework 
is remarkable mainly due to the possibility to ascribe stress-energy tensor even to 
low-frequency gravitational waves influencing the background, which is in contrast to 
standard linearization approach where the background is kept fixed. This is achieved 
by introducing a natural low-frequency cut-off, employing three different metrics (the 
prcmetric, the complete physical metric, and the average metric) and careful analysis 
of their mutual relations. 

On the other hand, in a now classic paper [3] Isaacson (inspired by previous 
works [4,5]) presented a perturbation method which can be used for studies of high- 
frequency gravitational waves. Such waves also influence the cosmological background 
in which they propagate. Isaacson's work stimulated further contributions in which 
his method was re- formulated using various formalisms, and explicitly applied to 
particular spacetimes, see e.g. [6-15]. 

In our present work we first briefly summarize and compare the two above 
mentioned perturbation schemes. In particular, it is shown that the Efroimsky 
method is not consistent if high-frequency gravitational waves are considered. Next 
(in section 3), we propose a possible modification of the Efroimsky formalism which 
may resolve this drawback. 

2. The formalism 

Efroimsky's approach [1, 2] is based on introducing three different smooth, non- 
degenerate, symmetric metrics on a differentiable manifold M, namely: 

(i) 7^^ — the "premetric" : vacuum metric corresponding to initial pure background 
without gravitational waves, 

(ii) <7 M „ — the "physical metric": full vacuum metric which describes both the 
background and the waves, 

(iii) q^ — the "average metric" : non- vacuum metric representing the background 
plus its perturbations with wavelength greater than L. In fact, it is the averaged 
full metric g llUl where the cut-off value L depends on observer's experimental 
abilities. Since no detector can measure gravitational waves of arbitrarily long 
wavelengths, the existence of such low-frequency cut-off is a natural assumption. 

One motivation for using these three distinct metrics is to resolve a (slight) discrepancy 
in standard linearization approach which considers only the metrics 7^, g^ v , and 
decomposition g^ v = 7 M , y + h^ v , where h^ v is a small perturbation. The contravariant 
components obtained as an inverse of g M „ are g^ v — -f^ — h^ v + 0(h 2 ), but 7^ is 
commonly used for raising and lowering indices. It is thus not clear which semi- 
Ricmannian manifold does this equality relate to. Such inconsistency can be ignored 
in the lowest order because it leads to the correct linear approximation of the wave 
equation. To extend the weak-field formalism to higher-order terms, the distinction 
between the premetric 7 and the average metric q is necessary, as it exhibits the back- 
reaction of the waves on the background geometry (Here and hereafter, indices of the 
metric tensors are sometimes suppressed for notational simplicity.) 
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Next step is to define the Ricci and Einstein tensors for an arbitrary metric g as 
-K-fj,i/{g) = [2 9 \9pv,[i + 3pM,f — 9fj.v,p)\,i ~ [29 [gp-y.p t 9pp,j ~ 9pi,p)\,v 

+ [^g lS (9ps n + g P1 ,s - g 1 s, P )\ [\g Sp {g P v^ + g Pl i,u - sw)] C 1 ) 

- [\g lp {g P s,v + g P u,s - g^s^)}[^g Sp (g P1 .,p, + g P p. n - g m , P )\ . 
Gpu{g) = Rpu{g) - \g^g al3 R a i3{g) , 

where g pT — (g)Zf, the same expressions apply to 7 and q. These equations remain a 
tensor even if we transfer to another semi-Riemann space (the reason is that covariant 
tensors are defined on a metric space rather than on some particular semi-Riemann 
one). From the proposals (i)-(iii) it follows that G M „(7) = = G l _ LU {g) , G^ v {q) ^ 0. 

Now, the differences between the covariant components of the above metrics are 
introduced, 

fipv = 9pv 9pv i \^) 

^]pv = qpu ^ipv • 

It is necessary to specify the semi-Riemann space: for raising or lowering indices 
and for covariant differentiation the averaged non-vacuum metric q will be used. 
Consequently, h and 77 are tensor fields on the semi-Riemann manifold (M, q), i.e. 

W = q pa q vP Kp , (3) 

if = g^j'V . 

Treating h^ as a perturbation of the metric q^ v the Ricci tensor Q) can be expanded 
in a power series 

Rpu (.9) = R$ (q) + R$ (q, h) + R$ (q, h) + R$ (q, h) + 0(h A ) , (4) 

where 

R^(q)=R flv (q), 

■E^ui/ \Qi "J = 2 |_2 \p,*^pT\v ~r ft {/^pr-.piy ~r flp,u\pr ^Tp\vp ^rv^p) \ 

h T ;p (h _ h \ _ (hpr -lhP-i r \(h -\-h —h )] 

lt> 1/ \' l Tp\p 'vpfJ.'jTJ \ llj \p 2 P J v b Tp\v ~ IL Tb>;fl " / fll/;Tjj • 

Analogously, 

^(7) = R$(q) + R$(q, (-V)) + R$(q, (-V)) + 0(rf) . (6) 

It is obvious that R^(q, (—rj)) = -R\j,v{q,rj) and R^(q, (—77)) = Rfiu(q,rf)- Accord- 
ing to assumptions that both g and 7 are vacuum metrics the following relation holds 

= #„„(«/) -iZ^ (7) 

= R$ (q, h) + R$ (q, h) + R$ (q, rj) + R$ (q, h) + 0{h A ) + 0( V 2 ) . (7) 

At this point Efroimsky sets three assumptions: 

Assumption 1. The perturbations h and rj are small in the sense that the terms 
of the orders 0(h 4 ) and 0(ry 2 ) are negligible. 

Assumption 2. The perturbations 77 and h 2 are of the same order. 

Assumption 3. The tensor field h consists of modes with short wavelengths which 
do not exceed the given maximal value L. 
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A physical interpretation of the perturbations given by is thus the following: hn V 
characterizes measurable gravitational waves whereas 7/ p „ is a shift of the background 
geometry from vacuum premetric 7 to nonvacuum effective average metric q due to 
the presence of gravitational waves. This enables us to interpret the equation Q as 
the wave equation for perturbations h on the background q — 7 + r\. To make this 
wave equation applicable, one has to express 77 in terms of h. Using the Brill-Hartle 
averaging procedure [5] over a spacetime volume of size L for J2J (Efroimsky considers 
only space averaging but when the measurement lasts much longer than the period of 
waves one can employ spacetime average) we obtain 

R$(q,V) = -(R$(q,h)) L . (8) 

The averaging brackets on the left-hand side are omitted because the term contains 
only the modes with wavelength greater than L. It is thus clear from JHJ and (0 
that the Assumption 2. is natural since the left-hand side is linear in 77 whereas the 
right-hand side is quadratic in h. 

Let us finally recall the derivation of the stress-energy tensor of gravitational 
waves. By analogy with the Ricci tensor expansion (0J the Einstein tensor of the 
vacuum premetric 7 is represented as a series 

= G^i) = G^{q) + G«(«, (-77)) + 0{if) , (9) 

and the effective stress-energy tensor of gravitational waves is defined as 

G^{q) = SttT^ = R$(Q,V) ~ hl^^R^V) • (10) 

From JHJl it follows (considering the Brill-Hartle averaging) that this tensor fully agrees 
with that of Isaacson [3] . 

The main advantage of the above Efroimsky's perturbation method is the 
possibility to consistently treat all low-frequency gravitational waves, and to explicitly 
derive effective stress-energy tensor (influencing the background) in this case. It can 
be extended to non-vacuum spacetimes with T M „ of ideal fluid and/or with a possible 
cosmological constant A, see [1,2]. However, there are some problems concerning 
high-frequency gravitational waves which will now be discussed. 

3. Modification to include high-frequency waves 

In this section we first explicitly demonstrate that one can not consistently 
apply Efroimsky's treatment on Isaacson's high-frequency waves [3] because the 
Assumption 2. is not fulfilled in such a case. Then we will present possible solution 
to this problem. 

Let us start with observation that it is the nonvacuum background curved by 
the presence of gravitational waves — not the vacuum premetric 7 — which is the 
basis of Isaacson's non-linear approach. Therefore, the nonvacuum average metric 
q is considered as the background on which high-frequency gravitational waves h 
propagate. 

We wish to use the Efroimsky formalism in the high-frequency regime such that 
the tensor field h contains high-frequency modes. We assume that they have short 
wavelengths A, and a small amplitude h = O(e), where e = A/5 < 1 is a small 
parameter because A <C S, S denoting a typical scale on which the background changes 
substantially. 

Let us emphasize that we follow here the same definition of the symbol 0(e n ) 
as in [3], namely / = 0(e n ) if there exists a constant C > such that |/| < Ce n 
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as £ — > 0. The quantity / need not necessarily be proportional to s n , it can be 
even smaller than Ce n for e — > 0. Therefore, the assumption /i = 0(e) does not 
automatically imply that h ~ e. The spectrum of possible high-frequency waves is 
thus not a priori restricted, it is only required that their amplitudes fall to zero at 
least linearly with e, i.e. \h(e)\ < Ce. 

Since we can consider S = O(l) it follows that O(e) = 0(X) anddh ~ h/X — 0(1). 
In accordance with Isaacson's approach (note that the decomposition now reads 
g = q + h, instead of the notation g = 7 + h used in [3]) we obtain the following 
orders of magnitude for the derivatives of the background q and the perturbation h: 

q^v = 0(1) , h^ v = 0(e) , 

q^a = 0(1) , V,« = 0(1) , (11) 

qnv, a p = O(l) , h^ap = 0(e _1 ) . 

This results in the orders of magnitude of the terms in the Ricci tensor expansion QJ j 
J5J as 

R$ = 0(l), R$ = 0(e~ 1 ), J$)=0(1), R$ = 0(e). (12) 

To apply the Efroimsky approach in this case we must consider the decomposition 
q = 7 + 77, where 7 is the vacuum premetric and 77 represents (in this case) substantial 
shift of the background geometry due to the presence of high-frequency gravitational 
waves h. We also introduce the scale L, such that J « i « S. This enables us 
simultaneously to consider an averaging procedure in accordance with the Isaacson 
approach, and also to introduce a meaningful cut-off scale L even if the wavelengths 
of high-frequency waves are not assumed to reach this value. 

Of course, the geometry shift r\ does not contain high-frequency perturbations. 
Considering the wave equation Q and using the Brill-Hartle averaging over a 
spacetime volume L to obtain the equation JSJ) we get in a conflict with the 
Assumption 2. which prescribes 0(n) — 0(h 2 ). Indeed, if h = 0(e) there should 
be 77 = 0(e 2 ). But the right-hand side of IJHJ) is now of the order O(l), see IJ12JI . and 
the same magnitude should also have the left-hand side. Since 77 does not contain 
high-frequency waves, it is essential that 77 = 0(1). This is obviously in contradiction 
with both Assumptions 1. and 2. In fact, it disables any consistent perturbation 
expansions in the powers of 77. 

Let us now suggest a modification of the Efroimsky formalism which will 
incorporate also the above case of a "substantial" change of the background geometry 
due to the presence of high-frequency waves. Instead of the perturbation expansion 
10 we consider a formal decomposition of the Ricci tensor of the premetric 7 = q — 77, 
namely 

= R^d) = R„M + A ^-(<7> (-»?)) > ( 13 ) 

by which equation the expression Ai? M „ is defined. Both terms on the right-hand side 
of 113fl are of the same order O(l). Moreover, the quantity AR^ is conserved with 
respect to the background geometry q which is easily seen from the equation 113|) and 
the relation (R fiu (q))' 1 ' — (the differentiation relates to the background metric q). 

The question concerning the gauge invariance of AR,^ with respect to generalized 
gauge transformations has been recently analyzed in detail by Anderson [16] in 
connection with possible definitions of the wave equation and stress-energy tensor 
for gravitational waves. Let us consider an arbitrary coordinate transformation of the 
type 

2" = a** + f" , (14) 



The Efroimsky formalism adapted to high-frequency perturbations 6 

that does not change the functional form of the background geometry q, i.e. 
q(x) — q(x) so that 7(2;) — > j{x) — q(x) — rj(x). Now, to prove the invariance of Ai? M „ 
we adopt (a slightly modified) Anderson's argumentation. Performing the above 
coordinate transformation (|14|l of the Ricci tensor decomposition (|13[1 we obtain 

R^(q(x)) + AR^(q(x), (-rim) = R^(l(x)) = . (15) 

Here R^ and AR^ are the same as R^ and AR^ U , respectively, because the 
definition (Q is maintained in any coordinates. Evaluating the relation (|15|l at x = x 
we thus get R llu (q(x)) = —AR fJ _^(q(x), (—rj(x))), and using l|13[) we obtain 

AR^(q(x), (-ri(x))) = AR^(q(x), (-rj(x))) . (16) 

A generalized gauge transformation is defined in [16] as a transformation in which 
the quantity rj(x) is substituted for n(x) into the tensor expressions of interest. This 
incorporates, as a particular case, the well-known infinitesimal gauge transformation, 

V„ v {%) = ?/> {x) + i^-.y + iv^ , (17) 

where 77, £ and their derivatives are small. Obviously, the equation l|16|l expresses a 
generalized gauge invariance of AR^ V . 

After introducing the above decomposition (|13fl and demonstrating its invariance 
we can now present modification and generalization of the Efroimsky formalism. 
Replacing the term i?}J(g,77) by —AR^ v (q,(—rj)) in equations JTJl, J8J, (fTTTfl . and 
omitting the terms 0(r/ 2 ) we obtain relations 

R$(q, h) + R${q, h) - AR pv (q, (-77)) + R$(q, h) + O^) = , (18) 

AR^(q,(- V )) = (R$(q,h)) L , (19) 

G^(q) = Sirf^ = -AR^(q, (-77)) + \q^ AR a0 (q, (-77)) . (20) 

In case when gravitational waves do not have high-frequency modes it is still possible 
to employ the expansion of —AR^ v (q, (—r))) in powers of r\ and use its dominant term 
RlJ(q,rj) instead. Thus we recover Efroimsky's previous results, cf. J7J), I©, (|TU1) . 

In general, however, expressing 77 in terms of h from the equation 1)19(1 becomes 
an extremely difficult task because it is no longer a linear equation for 77. To overcome 
this problem we can use the equation i|19[l and substitute for Ai? M „ into the remaining 
equations l|18|) and (|20[1 . We obtain the relations 

R$(q, h) + R$(Q, h) - (R$(q, h)) L + R$(q, h) + 0{h A ) = , (21) 

-G^q) = (R$(<1, h)) L - \ q , v q^(R { ^{ qi h)) L ee -SttT^ . (22) 

The equation l|22l) is obviously in perfect accordance with the Isaacson result [3] which 
represents the background response to the presence of high-frequency gravitational 
waves, using the Brill-Hartlc averaging to introduce the effective stress-energy tensor 
T^ H for high-frequency gravitational waves. The equation l(21|l is the wave equation 
for perturbations h on the average metric q. In the highest order of high-frequency 
approximation this clearly reduces to R^J = which also fully reproduces Isaacson's 
result. Additional terms in i121ll can be used for study of nonlinear effects on the wave 
propagation. 

Notice finally another interesting consequence of the equation 1(19J) and the gauge 
invariance IjRifl of Ai? M „ . This directly guarantees gauge invariance of the stress-energy 
tensor T^J 1 defined in (1221) (in the highest order). Proof of this property was presented 
already in the classic work [3], using however much more complicated method. 
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4. Concluding remarks 

In our contribution we have compared the Efroimsky [1,2] and the Isaacson [3] self- 
consistent perturbation schemes which describe propagation of weak gravitational 
waves on a cosmological background. In both these approaches the background is 
influenced by the waves, i.e. the non-linear effects are taken into account. The classical 
Isaacson method applies to high-frequency waves. On the other hand, the Efroimsky 
formalism is applicable to low-frequency gravitational waves but does not admit the 
high-frequency limit. We have suggested a modification of the Efroimsky formalism 
by employing the gauge-invariant decomposition l|13|) of the Ricci tensor, introduced 
recently by Anderson [16]. The resulting generalized system of equations i|18|) - (|20fl 
fully recovers the Efroimsky results in the absence of high-frequency modes, in the 
high-frequency limit it reproduces Isaacson's formulae. 

Although we have considered here for simplicity only vacuum metrics 7^ and 
g^u, possible generalization to non- vacuum spacetimes is straightforward. In fact, 
Efroimsky already generalized his formalism to spacetimes with ideal-fluid-like matter 
and a cosmological term [1,2], in case of the Isaacson high-frequency approach this 
has been done recently in [15]. 
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